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I. INTRODUCTION 

Two sets of exact solutions to Einstein's field equa- 
tion of General Relativity are of paradigmatic impor- 
tance: The first set describes the gravitational field of 
quasi-isolated objects in an asymptotically flat space- 
time. Among them is the exterior Schwarzschild solution 
that describes the stationary gravitational field outside 
a spherically symmetric star or black hole of mass m 
with vanishing intrinsic angular momentum (spin) and 
vanishing electric charge. (The latter two features being 
included in the three-parameter Kerr-Newman family of 
solutions.) Such asymptotically flat solutions are meant 
to apply to a region outside the central object which, 
on one hand, must be sufficiently far from the considered 
object, so as to legitimately neglect small irregularities of 
its surface and/or small deviations from perfect spheri- 
cal symmetry. On the other hand, and more importantly, 
the region of applicability must also be sufficiently close 
to the considered object in order not to include, or come 
close to, other compact sources, or not contain too much 
dust-filled space between it and the object which would 
also act as disturbing source for the gravitational field. 
In particular, the large-distance asymptotic behavior of 
such solutions is an idealization and not meant to be 
strictly that of any object in the real world. 

On the other hand, the second set of paradig- 
matic solutions are the cosmological ones, which aim 
to model the behavior of space-time at the largest 
cosmological scales, without trying to be realistic at 
smaller scales. Among them is the family of homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) cosmologies on which the cosmological 
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standard-model is based. 

Given that situation, the task is to combine the virtues 
of both classes of solutions without the corresponding 
deficiencies. This means to find exact solutions for the 
gravitational field of a compact object 'immersed' (see 
below) into an otherwise cosmological background. This 
would appear to be an easy task if the field equations 
were linear, for, in that case, one would just add the so- 
lution that describes the gravitational field of a compact 
object in an otherwise empty universe to the cosmological 
solution that corresponds to a homogeneous distribution 
of background matter. Here the mathematical opera- 
tion of addition appears to be the obvious realization of 
what one might be tempted to call 'simultaneous physical 
presence' and hence, in view of the individual interpre- 
tations of both solutions, the 'immersion' of the compact 
object into the cosmological background. But this imme- 
diate interpretation in physical terms of a simple mathe- 
matical operation is deceptive. This becomes obvious in 
non-linear theories, like General Relativity (GR), where 
no simple mathematical operation exists that produces 
a new solution out of two old ones and where the very 
same physical question may still be asked. 

The proper requirement for a mathematical represen- 
tation of the envisaged physical situation must, first of 
all, consist in asymptotic conditions which ensure that 
the sought for solution approximates the given (e.g. 
Schwarzschild) one for small distances and a particu- 
lar cosmological one (e.g. FLRW) for large distances. 
Second, it must specify somehow the physics in the in- 
termediate region. Usually this will include a specifica- 
tion of the matter components and their dynamical laws 
together with certain initial and boundary conditions. 
Needless to say that this will generally result in a complex 
system of partial differential equations. Most analytic 
approaches therefore impose further simplifying assump- 
tions that automatically guarantee the right asymptotic 
behavior and at the same time reduce the free functions 
to a manageable number. 

In this paper we will discuss a particular such ap- 
proach, which is originally due to McVittie [l^ and which 
has been further analyzed and clarified in a series of care- 
fully written papers by Nolan [l^, [l^, Our main 
motivation is that recently McVittie's solution has been 
severely criticized as not being able at all to model the 
envisaged situation[3,Q, whereas a family of slightly gen- 
eralized ones in which some restrictions concerning 
the motion of matter and the existence of heat flows is 
lifted, is argued to be free of the alleged problems. The 
existence of an exact solution to Einstein's equation that 
models local inhomogeneities is clearly of great impor- 
tance, for example in estimating reliable upper bounds 
to the possible influence of global cosmological expansion 
onto the dynamics and kinematics of local systems jlj . 

The paper is organized as follows: In SectionlU we re- 
view what we call the McVittie model. We discuss its 
metric ansatz and what its entails regarding the geom- 
etry of spacetime. Then we discuss the assumptions re- 



garding the motion of the matter and how this, together 
with Einstein's equation, determines one of the two free 
functions in the metric ansatz as a simple function of the 
other. We interpret this condition in terms of an appro- 
priate concept of local mass as saying that the object 
does not accrete mass from the ambient matter. In Sec- 
tion lllll we take a second and closer look at the McVittie 
ansatz and note some of its characteristic features which, 
we feel, have not sufficiently carefully been taken into ac- 
count in d, 0, 0] . In the light of these observations we 
then discuss in Section lTVl the attempted generalizations 
of McVittie's solution in the references just mentioned. 
We find that some of the conclusions drawn are indeed 
unwarranted. 



II. THE MCVITTIE MODEL 

The characterization of the McVittie model is made 
through two sets of a priori specifications. The first set 
concerns the metric (left side of Einstein's equations) and 
the second set the matter (right side of Einstein's equa- 
tions). The former consists in an ansatz for the metric, 
which can formally be described as follows: Write down 
the Schwarzschild metric for the mass parameter m in 
isotropic coordinates, add a conformal factor a^{t) to the 
spatial part, and allow the mass parameter m to depend 
on time. Hence the metric reads 



I - m{t)/2r 
l + m{t)/2r 

4 



m 



it) 



2r 



dt^ 



\t) (dr^+r^g,.).. 



(1) 



where g<,2 = d9^ + sin^ 9dip^ is the standard metric on 
the unit 2-sphere. Here we restricted attention to the 
asymptotically spatially flat (i.e. k — 0) FLRW metric, 
which is compatible with current cosmological data [ll|. 
For simplicity we shall refer to (U) simply as McVittie's 
ansatz, though this is not quite correct since McVittie 
started from a general spherically symmetric form and 
arrived at (H]) with m{t)a{t) — const, after imposing a 
condition that he interpreted as condition for no matter- 
infall. The ansatz ([T]) is obviously spherically symmetric 
with the spheres of constant radius r being the orbits of 
the rotation group. "'^ In the next section we will discuss 
in more detail the geometric implications of this ansatz, 
independent of whether Einstein's equation holds. 



^ "Spherical symmetry" of a spacetime means the following: There 
exists an action of the group 50(3) on spacetime by isometrics, 
which is such that the orbits are either two-dimensional and 
spacelike or fixed points. The "spheres" implicitly referred to in 
this term correspond to the two-dimensional orbits, even though 
they might in principle also be two-dimensional real projective 
spaces. In the cases we discuss here they will be 2-spheres. 
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As already discussed in the introduction, the model 
here is meant to interpolate between the spherically sym- 
metric gravitational field of a compact object and the 
environment. It is not to be taken too seriously in the 
region very close to the central object, where the basic 
assumptions on the behavior of matter definitely turn 
unphysical. However, as discussed in at radii much 
larger than (in geometric units) the central mass (to be 
defined below) the k = McVittie solution seems to pro- 
vide a viable approximation for the envisaged situation. 

The second set of specifications, concerning the matter, 
is as follows: The matter is a perfect fluid with density 
g and isotropic pressure p. Hence its energy-momentum 
tensor is given by^ 



gM(8)U+p(M(8)tt — g) 



(2) 



Furthermore, and this is where the two sets of specifi- 
cations make contact, the motion of the matter (i.e. its 
four-velocity field) is given by 



u ^ Bq . 



(3) 



where eg is the normalization of d/dt (compare (|10l) '). 
Finally, the explicit cosmological constant on the left- 
hand side of Einstein's equation is assumed to be zero, 
which implies no loss of generality, since a non-zero cos- 
mological constant can always be regarded as special part 
of the matter's energy-momentum tensor (compare lIVp . 
No further assumptions are made. In particular, an equa- 
tion of state, like p = p{g), is not assumed. The reason 
for this will become clear soon. Later generalizations will 
mainly concern ([2]) and ([3]). 

The Einstein equation-^ now links the specifications 
of geometry with that of matter. It is equivalent to 
the following three relations between the four functions 
m{t), a(t), g{t, r), and p(t, r): 



{a m) 
Sng 

Sirp 



0, 









) 








d 


-( 






a 



1 + m/2r 
1 — m/2r 



(4a) 
(4b) 

(4c) 



Note that here Einstein's equation has only three inde- 
pendent components (as opposed to four for a general 
spherically symmetric metric), which is a consequence 
of the fact that the Einstein tensor for the McVittie 
ansatz ([ij is spatially isotropic. This will be discussed in 
more detail in the next section. 



^ Here and in what follows we denote the metric-dual (1-form) 
of a vector u by underlining it, that is, u := g{u, ■ ) is the 1- 
form metric-dual to the vector u. In local coordinates we have 
u = Si'' 9^ and u = Mpdx^', where ti^ := g^^u". 

^ We speak of "the Einstein equation" in the singular since we 
think of it as a single tensor equation, which only upon introduc- 
ing a coordinate system decomposes in many scalar equations. 



Equation (|4ap can be immediately integrated: 
TO(i) 



(5) 



where mo is an integration constant. Below we will show 
that this integration constant is to be interpreted as the 
mass of the central body. 

Clearly the system Q is under-determining. This is 
expected since no equation of state has yet been imposed. 
The reason why we did not impose such a condition can 
now be easily inferred from ([4]): whereas (j4bp implies 
that g only depends on t, (f4c|) implies that p depends 
on t and r iff {a/a)' ^ 0. Hence a non-trivial relation 
p — p{g) is simply incompatible with the assumptions 
made so far. The only possible ways to specify p are 
p = or + p = 0. In the first case (Pc| implies that 
aja = if 7no 7^ (since then the second term on the 
right-hand side is r dependent, whereas the first is not, 
so that both must vanish separately), which corresponds 
to the exterior Schwarzschild solution, or a{t) oc t^^^ if 
mo — 0, which leads to the fiat FLRW solution with 
dust. In the second case the fluid just acts like a cos- 
mological constant A — Sng (using the equation of state 
g + p — in div T = it implies dp = and this, in 
turn, using again the equation of state, implies dg — 0) 
so that this case reduces to the Schwarzschild-de Sitter 
solution. To see this explicitly, notice first that (|4bl4cp 
imply the constancy of -ff = a/ a = \J A/3 and hence 
one has a{t) — oq exp (t y^A/s) . With such a scale- 
factor the McVittie metric ^ with ([5]) turns into the 
Schwarzschild-de Sitter metric in spatially isotropic coor- 
dinates. The explicit formulae for the coordinate trans- 
formation which brings the latter in the familiar form can 
be found in Section 5 of [3 and also in Section 7 of • 
Finally, note from (Pa|) that constancy of one of the func- 
tions m and a implies constancy of the other. In this case 
()4bl4cp imply p = p = 0, so that we are dealing with the 
exterior Schwarzschild spacetime. 

A specific McVittie solution can be obtained by choos- 
ing a function a(t), corresponding to the scale function 
of the FLRW spacetime which the McVittie model is 
required to approach at spatial infinity, and the con- 
stant mo, corresponding to the 'central mass'. Rela- 
tions (|4bl4cp . and ([5]) are then used to determine g, p, 
and m, respectively. Clearly this 'poor man's way' to 
solve Einstein's equation holds the danger of arriving 
at unrealistic spacetime dependent relations between g 
and p. This must be kept in mind when proceeding in 
this fashion. For further discussion of this point we refer 

to[il,p. 

As will be discussed in more detail in Section llll CI be- 
low, in the spherically-symmetric case the concept of lo- 
cal mass (or energy) is well captured by the Misner-Sharp 
(MS) energy [l3|, whose purely geometric definition in 
terms of Riemannian curvature allows to decompose it 
into a sum of two terms, one of which comes from the 
Ricci- the other from the Weyl curvature. It is the latter 
which may be identified with the gravitational mass of 
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the central object. Applied to ([T]), the Weyl contribution 
to the MS energy can be written in the following form, 
also taking into account 

47r 

£^R=Y^?'^^, (6a) 

i?w — m-o . (6b) 

The constancy of i?w is then interpreted as saying that 
no energy is accreted from the ambient matter onto the 
central object. 

We now briefly discuss the basic properties of the mo- 
tion of cosmological matter. Being spherically symmet- 
ric, the velocity field u specified in ^ is automatically 
vorticity free. The last property is manifest from its hy- 
persurface orthogonality, which is immediate from ([T]). 
Moreover, u is also shear free. This, too, can be imme- 
diately read off ^ once one takes into account the fol- 
lowing result, whose proof we sketch in Appendix[Cl A 
spherically symmetric normalized timelike vector field u 
in a spherically symmetric spacetime {M,g) is shear free 
iff its corresponding spatial metric, that is, the metric g 
restricted to the subbundle := {v £ TA4 \ g{v, u) = 
0}, is conformally flat. The metric ([T]) obviously is spa- 
tially conformally flat with respect to the choice ^ made 
here. Moreover, the expansion (divergence) of u is 

e = 3H, (7) 

where H := a /a, just as in the FLRW case. In par- 
ticular, the expansion of the cosmological fluid is homo- 
geneous in space. Exactly as in the FLRW case is also 
the expression for the variation of the areal radius along 
the integral lines of u (that is the velocity of cosmologi- 
cal matter measured in terms of its proper time and the 
areal radius): 

u{R) = HR, (8) 

which is nothing but Hubble's law. Recall that for a 
spherically symmetric spacetime the areal radius, de- 
noted here by R, is the function defined by R{p) := 
y^4(p)/47r, where A{p) is the proper area of the 2- 
dimensional S'0(3)-orbit through the point p. For the 
McVittie spacetime the areal radius is given explicitly 
in The acceleration of u, which in contrast to the 

FLRW case does not vanish here, is given by 

_ Too fl + m/2r\ 
^-^ = l^[l^r)^^- 

Here ei is the normalized vector field in radial direction 
as defined in (jlOp . In leading order in uiq/R this corre- 
sponds to the acceleration of the observers moving along 
the timelike Killing field in Schwarzschild spacetime. 

It is also important to note that the central gravita- 
tional mass in McVittie's spacetime may be modeled by a 
shear-free perfect-fluid star of positive homogeneous en- 
ergy density p^ . The matching is performed along a 
world-tube comoving with the cosmological fluid, across 



which the energy density jumps discontinuously. This 
means that the star's surface is comoving with the cos- 
mological fluid and hence, in view of (O, that it geomet- 
rically expands (or contracts). This feature, however, 
should be merely seen as an artifact of the McVittie 
model (in which the relation ([7]) holds), rather than a 
general property of compact objects in any cosmological 
spacetimes. Positive pressure within the star seems to 
be only possible if 2ad -I- < (see Eq. (3.27) in [l4| 
with a = exp(/3/2)), that is, for deceleration parameters 
g> 1/2. 

III. GEOMETRY OF THE MCVITTIE ANSATZ 

In this section we will discuss the geometry of the met- 
ric IH) independent of the later restriction that it will 
have to satisfy Einstein's equation for some reasonable 
energy-momentum tensor. This means that at this point 
we shall not assume any relation between the two func- 
tions m{t) and a{t), apart from the first being non neg- 
ative and the second being strictly positive. We will 
discuss the metric's 'spatial Ricci-isotropy' (a term ex- 
plained below), its singularities and trapped regions, and 
also compute its Misner-Sharp energy decomposed into 
the Ricci and Weyl parts. We shall start, however, by an- 
swering the question of what the overlap is between the 
geometries represented by ((J) and the conformal equiva- 
lence class of the exterior Schwarzschild geometry. 

A. Relation to conformal Schwarzschild class 

This question is an obvious one in view of the way 
in which ([T]) is obtained from the exterior Schwarzschild 
metric. It is clear that for to = toq = const, the 
metric (HI) is conformally equivalent to the exterior 
Schwarzschild metric, since upon using a new time co- 
ordinate T with dT = dt/a(t) we can pull out a^(t) as 
a common conformal factor. The following proposition, 
whose proof wc shall give in Appendix [XI states that a 
constant to is in fact also necessary condition: 

Proposition 1. Let iSmcv denote the set of metrics in 
the form of the McVittie ansatz {Ip (parametrized by the 
two positive functions a and to ) and S^s the set of met- 
rics conformally equivalent to an exterior Schwarzschild 
metric (parametrized by a positive conformal factor and 
a constant positive Schwarzschild mass Mq). Then the 
intersection between S^cv eind S^s given by the subset 
of metrics in iSmcv with constant to or, equivalently, by 
the subset of metrics in S^s whose conformal factor has 
a gradient proportional to the Killing field d / dT of the 
Schwarzschild metric (see iAlb\) for notation). 

Note that we excluded the 'trivial' cases in which m or 
Mq (or both) vanish for the following reason: Comparing 
the expressions for the Weyl part of the MS energy of 
the two types of metrics (see (jA10[) in Appendix [X| it 
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follows that m vanishes iff Mq does and this, in turn, 
leads to a metric conformally related to the Minkowski 
metric where the conformal factor depends only on time, 
that is, a FLRW metric. But such a spacetime, being 
homogeneous, is not of interest to us here. 

In particular, Proposition [T] implies that the metric of 
Sultana and Dyer 19,] are not of type as suggested 
in Section IV A of ^Jj and allegedly shown in Section II 
of 0] (cf. our footnote [S] at page [5]). This immediately 
follows from the observation that the conformal factor, 
expressed as function of the standard Schwarzschild co- 
ordinates that appear in (jAlb[) . is given by ri(r, i?) — 
(T + 2Moln(i?/2Mo - 1))^ (compare Eqs. (8) and (9) 
of [i^), which also depends on R and hence does not 
satisfy the condition of Proposition [TJ We will have to 
say more about this at the beginning of Section lTvl and 
in SectionHVD] 



B. Spatial Ricci-isotropy 

An important feature of any metric that is covered 
by the ansatz ([T]) is, that its Einstein tensor is spatially 
isotropic in the following sense: 'Spatially' refers to the 
directions orthogonal to d/dt and 'isotropy' to the con- 
dition that the spatial restriction of the spacetime's Ein- 
stein tensor is proportional to the spatial restriction of 
the metric. Note that, since the spacetime's metric is 
time dependent, the spatial restriction of the spacetime's 
Einstein or Ricci tensor is not the same as the Einstein 
or Ricci tensor of the spatial sections with their induced 
metrics. Hence the notion of spatial isotropy of the Ein- 
stein tensor used here is not the same as saying that the 
induced metric of the slices is an Einstein metric. 

Given that the Einstein tensor of ([T]) is spatially 
isotropic in the sense used here, it is then obvious that 
Einstein's equation will impose a severe restriction upon 
the matter's energy-momentum tensor, saying that it, 
too, must be spatially isotropic. The degree of special- 
ization implied by this will be discussed in more detail 
below. Here we only remark that this observation al- 
ready answers in the negative a question addressed, and 
left open, in the last paragraph of [1], of whether ([1]) 
is the most general spherically symmetric solution de- 
scribing a black hole embedded in a spatially flat FLRW 
background: It clearly is not. 

In passing we make the obvious remark that, since 
Ein = Ric — (1/2) Scalg, where Ric denotes the Ricci 
tensor, the Einstein tensor is spatially isotropic iff the 
same holds for the Ricci tensor. For this reason we will 
from now on refer to spatial Ricci-isotropy to denote the 
feature in question. 

Now, a way to actually show spatial Ricci-isotropy is 
to compute the components of the Einstein tensor with 
respect to the orthonormal tetrad {e^lpg^o,- - ,3} of ([T]) 
defined by 



where {x^} = {t, r, 9, <p}. Here, and henceforth, we write 
\\v\\ ^y\g{v, v)\. Note that eo,ei are orthogonal to 
and 62, 63 tangent to the 2-spheres of constant radius r. 
The non-vanishing independent components of the Ein- 
stein tensor with respect to the orthonormal basis (fTT 



Ein(eo,eo) 
Ein(eo,ei) 
Ein(ej, Sj) 



= if(4)'(«"^)'' 

+ 24f] Si 



(11a) 
(lib) 
(11c) 



where an overdot denotes differentiation along d/dt. Be- 
fore explaining the functions A, B, R, and F, note that 
the spatial isotropy of the Einstein tensor follows imme- 
diately from (|llcp . since Ein(ei,ej) cx Sij. In pT|) and 
in the following we set: 

A{t,r) ■.= l + m{t)/2r, B{t,r) := 1 - m{t) /2r , (12) 

and 



R{t,r)=[l + ^) a{t)r. 



(13) 



where R is the areal radius for the McVittie ansatz ([T]), 
and also 



1 (am)' 
rB a 



(14) 



In passing we note that both quantities, F and a m, that 
appear in the components of the Einstein tensor, have 
a geometrical interpretation: the former is one third the 
expansion of the vector field Cq, that is, F = div(eo)/3, 
and the latter is the Weyl part of the Misner-Sharp en- 
ergy of the metric ([1]) (see (f22|) . below). Moreover, as we 
already noted in Section |TT1 the observer field Cq is free of 
vorticity and shear. Hence, taking into account the rela- 
tion (jCSP between the expansion 6 and the shear scalar a 
of an arbitrary spherically-symmetric observer field, the 
expansion of Cq can be simply written as 3dR{eQ) / R so 
that F may be expressed as 



F = dR{eo)/R. 



(15) 



In order to estimate the degree of specialization im- 
plied by spatial Ricci-isotropy, we ask for the most gen- 
eral spherically symmetric metric for which this is the 
case. To answer this, we first note that any spherically 
symmetric metric can always be written in the form 



B{t,r) 



dt^ 



(10) 



^^^^^^ , a^{t)A\t.r){dr^+r'gs^). (16) 

This reduces to McVittie's ansatz ^ \i A,B are given 
by (|12p. For the general spherically symmetric metric 
(I16p . spatial Ricci-isotropy can be shown to be equivalent 
to 

5'^{AB)-S{5A){6B)^Q, (17) 

where 5 := r^^d/dr = 2d/dr^. It is obvious that there 
are many more solutions to this differential equation than 
just (HH). 
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C. Misner— Sharp energy 

In order to be able to interpret ([T]) as an ansatz for an 
inhomogeneity in a FLRW universe, it is useful to com- 
pute the Misner-Sharp (MS) energy and, in particular, 
its Ricci and Weyl parts. This concept of quasi-local 
mass, which is defined only for spherically symmetric 
spacetimes, and which in this case coincides with Hawk- 
ing's more general definition Q of quasi-local mass (see 
e-g- [H)i allows to detect localized sources of gravity. 

We recall the geometric definition of the MS energy [1, 



E := -^R^K, 



(18) 



where R denotes the areal radius and K the extrinsic 
curvature. More precisely, the equation should be read 
and understood as follows: First of all, the quantities 
R and K, and hence also E, are real- valued functions 
on spacetime. In order to determine their values at a 
point p, recall that, due to the requirement of spherical 
symmetry, there is a unique two- (or zero-) dimensional 
50(3) orbit S{p) through p. The value of _R at p is as 
explained below Eq. ([8]) and the value oi K at p is 



K{p) :-- 



Rieui{Xp,Yp,Xp,Yp) 



giXp, Xp)g{Yp, Yp) - {g{Xp,Yp)y 



(19) 



Here Riem is the (totally covariant) Riemannian cur- 
vature tensor of spacetime and Xp and Yp are any two 
linearly independent vectors in the tangent space at p 
which are also tangent to the orbit S{p). Note that the 
right-hand side only depends of the plane spanned by 
Xp, Yp and not on the vectors spanning it. Finally we 
note that the minus sign in is just a relict of our 
signature choice (mostly minus). 

From the curvature decomposition for a spherically 
symmetric metric (see [l|) one can rewrite (jl8p in the 
form 



R 



E=-{l+g{VR, VR)) 



(20) 



where Vi? denotes the gradient vector-field of R. This 
provides a convenient expression for the computation of 
the MS energy. For a self-contained review of the basic 
properties of the MS energy as well as its interpretation 
as the amount of active gravitational energy contained 
in the interior of the spheres of symmetry (50(3)-orbits) 
and its relation with the other mass concepts, see |l|- 

The decomposition of the Riemann tensor into a Ricci 
and a Weyl part leads, together with (|18p. to a natural 
decomposition of the MS energy into a Ricci and Weyl 
part (see also [l|). For the Ricci part of the MS energy 
of HI) we get 

E^ = Ei^^eo, eo) = ^{dR{eo)f . (21) 

The first equality in (|2ip can be derived by merely using 
the spatial Ricci-isotropy in the expression for the Ricci 



part of the Riemann tensor. The second equality follows 
then with (jllap and (HSj). The Weyl part can now be 
obtained as the difference between the full MS energy and 
PTjl . We use the expression ([^0]) for the former and write 

g{VR,VR) = {eo{R)y - {ei{R)f . The part involving 
Co {R) equals the Ricci part (pij) and hence the Weyl part 
is given by (i?/2)(l - (ei(i?))2). From (|13|) we calculate 
ei(i?) and hence obtain for the Weyl part of the MS 
energy: 



Em 



a m . 



(22) 



The Ricci part of the MS energy is that part which, via 
Einstein's equation, can be locally related to the mat- 
ter's energy-momentum tensor, whereas the gravitational 
mass of the central object is contained in the Weyl part 
of the MS energy. Notice that the latter is spatially con- 
stant (the functions a and m in (j22p only depend on time) 
but may depend on time. If the latter is the case we in- 
terpret this as saying that the central mass exchanges 
energy with the ambient matter. 



D. Singularities and trapped surfaces 

Next we comment on the singularity properties of the 
McVittie ansatz (|T]). From (|llcp one suspect, because of 
the term proportional to 1/B, a singularity in the Ricci 
part of the curvature at r = m/2 (that is at i? = 2am = 
2E^). In fact, this corresponds to a genuine curvature 
singularity, as one can see from looking, for example, at 
the following expression for the scalar curvature (i.e. the 
Ricci scalar). 



Seal: 



-12F^ 



6-F 



(23) 



which can be quickly computed from (jlip . In Appendix [Bl 
we insert into this expression the definition (jl4p of F and 
expand this in powers of 1/ {rB) . This allows to prove 

Proposition 2. The Ricci scalar for a metric of the form 
(Q]) becomes singular in the limit r ^ m/2 for any func- 
tions a and m, except for the following three special cases: 
(i) m — and a arbitrary (FLRW), 
(ii) a and m are constant (Schwarzschild) , and 
{Hi) {am)' = and {a/ a)' = (Schwarzschild~de Sitter). 

This means that, as long as we stick to the ansatz ([T|), at 
r — m/2 there will always (with the only exceptions listed 
above) be a singularity in the Ricci part of the curvature 
and thus, assuming Einstein's equation is satisfied, also in 
the energy momentum tensor, irrespectively of the details 
of the underlying matter model. Hence any attempt to 
eliminate this singularity by maintaining the ansatz ([1]) 
and merely modifying the matter model is doomed to 
fail. 

In particular, this is true for the generalizations pre- 
sented in Q, contrary to what is claimed in that work 
and its follow ups [1, @ . We also remark that it makes no 
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sense to absorb the singular factors 1 /B in front of the 
time derivatives by writing {A/ B)d / dt as eg and then 
argue, as was done in [3], that this eliminates the singu- 
larity. The point is simply that then cq applied to any 
continuously differentiablc function diverges as r — > m/2. 
Below we will argue that this singularity lies within a 
trapped region. 

Specializing to the McVittie model, recall that in this 
case it is assumed that the fluid moves along the integral 
curves of d/dt, which become lightlike in the limit as r 
tends to m/2. Their acceleration is given by the gradient 
of the pressure, which necessarily diverges in the limit 
r to/2, as one explicitly sees from For a more 
detailed study of the geometric singularity at r = m/2, 
see [Sill. 

For spherically symmetric spacetimes the Weyl part of 
the curvature has only a single independent component, 
which is —2/B? times the Weyl part of the MS energy, 
by the very definition of the latter (see [l|). The square 
of the Weyl tensor for the ansatz ([1]) may then be conve- 
niently expressed as 



by Rs can be written in the form 



(Weyl, Weyl) = 48 



(a m)^ 



(24) 



This shows that R — also corresponds to a genuine 
curvature singularity, though this is not part of the region 
covered by our original coordinate system, for which r > 
to/2 (that is R> 2E„). 

It is instructive to also determine the trapped regions 
of McVittie spacetime. Recall that a spacelike 2-sphere 
S is said to be trapped, marginally trapped, or untrapped 
if the product 9'^9~ of the expansions (for the defini- 
tions see e.g. [l|) for the ingoing and outgoing future- 
pointing null vector fields normal to S is positive, zero, 
or negative, respectively. Taking S to be Sr, that is, 
an 50(3) orbit with areal radius R, it immediately fol- 
lows from the relation 2 9+9- = g(Vi?, VR)/R^ (see Q) 
that Sr is trapped, marginally trapped, or untrapped iff 
g(Vi?, Vi?) is positive, zero, or negative, respectively. 
This corresponds to timelike, lightlike, or spacelike dR, 
or equivalently, in view of ((20)) . to 2E — R being posi- 
tive, zero, or negative, respectively. Using (PT|) together 
with piap . the MS energy for the McVittie ansatz can 
be written a,s E = E„ + i?^F^/2, so that 



2E-R = F^R^ -R + Rs. 



(25) 



Here we defined the 'Schwarzschild radius' as Rs := 2i?w, 
which generally will depend on time. We wish to de- 
termine the values of the radial coordinate (r or R) at 
which the expression ([^5]) assumes the value zero. We 
shall continue to work with R rather than r since R has 
the proper geometric meaning of areal radius. In the 
region we are considering (that is r > m/2 or, equiv- 
alently, R > Rs) the inversion of (fT3|) reads r{R) = 
R{1 - Rs/2R + ^\ - Rs/R)l2a, so that ^ divided 



2E-R 



Rs 



= \ T1 



- 1 + ^Jx{x - 1) 



x^-x + l. (26) 



Here we introduced the dimensionless radial coordinate 
X :— R/Rs and the (small) parameters e :~ Rs and 
T] :— Rs/Rh, where Rh '■= ^/H denotes the 'Hubble 
radius'. Recall that since R > Rs we have a; > 1. 

Consider first the McVittie case, in which e = Rs = 0. 
Then (j26p turns into a cubic polynomial in x which is 
positive for x = and tends to ±oo for x — s- ±oo. Hence 
it always has a negative zero (which does not interest us) 
and two positive zeros iff 



Rs/Rh < 2/3V3 « 0.38. 



(27) 



This clearly corresponds to the physical relevant case 
where the Schwarzschild radius is much smaller than 
the Hubble radius. One zero lies in the vicinity of 
the Schwarzschild radius and one in the vicinity of the 
Hubble radius, corresponding to two marginally trapped 
spheres. The exact expressions for the zeros can be easily 
written down, but are not very illuminating. In leading 
order in the small parameter rj — Rs/Rh, they are ap- 
proximated by 



i?i =i?s(l + r;2 + 0(r;4)), 
R2=RH{l-ri/2 + 0{ri^)) 



(28a) 
(28b) 



From this one sees that for the McVittie ansatz the radius 
of the marginally trapped sphere of Schwarzschild space- 
time {Rs) increases and that of the FLRW spacetime 
{Rh) decreases. The first feature can, for the McVit- 
tie model, be understood as an effect of the cosmological 
environment, whereas the latter is an effect of the inho- 
mogeneity in form of a central mass abundance. All the 
spheres with R < Ri or R > R2 are trapped and those 
with Ri < R < R2 are untrapped. In particular, the 
singularity r — m/2, that is i? = 2E^ — Rs, lies within 
the inner trapped region. 

In the case in which e = Rs is non-zero and 'small' 
(see below in which sense) , we expect that the zeros (pS)) 
vary smoothly in e so that, in particular, the singularity 
at i? = Rs still remains within the inner trapped region. 
An expansion in e gives, for the zero in the vicinity of 
the Schwarzschild radius: 

Ri (e) - i?s (1 + + (2 - 2?7 + 13r;2)e + 0{r^^s^)) , (29) 



which clearly reduces to (|28ap for e = 0. From this ex- 
pression one sees that, according to the physical expec- 
tation, in case of accretion {e > 0) the inner marginally 
trapped sphere becomes larger in area, whereas in the op- 
posite case (e < 0) it shrinks. In our approximation (|29p . 
the singularity R = Rs continues to lie inside the trapped 
region for 'accretion rates' e = Rs > — ry'^/2 or, in terms 
of physical quantities and re-introducing the factors of c. 
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for Rs/c > —{Rs/Rh)'^/'^- However, this also character- 
izes the region of validity of the expansion ([^5]) : Given 
a positive r/, an expansion in e around zero exists only 
for e > since there exists no expansion on both 

(e, rf) around (0, 0) (this is because the partial derivative 
of ([26]) with respect to x does not exist at a; = 1). 

E. Other global aspects 

Another aspect concerns the global behavior of the 
McVittie ansatz ([T]). We note that each hypersurface 
of constant time i is a complete Riemannian manifold, 
which, besides the rotational symmetry, admits a dis- 
crete isometry given in (r, 6, ip) coordinates by 

^{r,0,ip)^{im/2)^r-\d,ip). (30) 

This corresponds to an inversion at the 2-sphere r = 
to/2, which shows that the hypersurfaces of constant t 
can be thought of as two isometric asymptotically-flat 
pieces joined together at the 2-sphere r — m/2. This 2- 
sphere is totally geodesic since it is a fixed-point set of an 
isometry; in particular, it is a minimal surface. Except 
for the time-dependent factor m{t), this is just like for 
the slices of constant Killing time in the Schwarzschild 
metric (the difference being that ([30| does not extend to 
an isometry of the spacetime metric unless to = 0). Now, 
the fact that r ^ corresponds to an asymptotically flat 
end of each of the 3-manifolds t — const, implies that 
the McVittie metric cannot literally be interpreted as 
corresponding to a point particle sitting at r = (r = 
is in infinite metric distance) in an otherwise spatially 
flat FLRW universe, just like the Schwarzschild metric 
does not correspond to a point particle sitting at ?■ = 
in Minkowski space. Unfortunately, McVittie seems to 
have interpreted his solution in this fashion jl2| which 
even until recently gave rise to some confusion in the 
literature (e.g. [B,[3,[23|)- A clarification was given in p^ . 

IV. ATTEMPTS TO GENERALIZE MCVITTIE'S 
MODEL 

The first obvious generalization consists in allowing for 
a non-vanishing cosmological constant. However, as was 
already indicated before, this is rather trivial since it 
merely corresponds to the substitutions p — > g + qa and 
p ^ p + PA in jll, where qa ■— A/Stt and pA ■= — A/Stt 
are the energy-density and pressure associated to the cos- 
mological constant A. 

The attempts to non-trivially generalize the McVittie 
solution have focused so far on keeping the ansatz ^ and 
relaxing the conditions on the matter in various ways. In 
Q generalization were presented allowing radial fiuid mo- 
tions relative to the observer vector field d/dt (that is 
relaxing condition as well as including heat conduc- 
tion. Below we will critically review these attempts, tak- 
ing due care of the geometric constraints imposed by the 



ansatz ([T]), and also outline how to explicitly construct 
the respective solutions. 



Another exact solution that models an inhomogene- 
ity in a cosmological spacetime was presented in [l9| 
by Sultana and Dyer and was recently analyzed in 
Here the metric is conformally equivalent to the exte- 
rior Schwarzschild metric and the cosmological matter is 
composed of two non-interacting perfect fluids, one being 
pressureless dust, the other being a null fluid. One might 
ask if this solution flts into the class of McVittie models, 
as was suggested in |4j* and allegedly conflrmed explicitly 
in f^l^. However, as we already noted at the end of Sec- 
tion lHI Al above in view of Proposition[l] this is not the 
case. Two further way to see this are as follows: First, 
the Sultana-Dyer metric is not spatially Ricci-isotropic^ 
and, second, the McVittie metric is not compatible with 
the matter model used by Sultana and Dyer, with the 
sole exception of trivial or exotic cases, as will be shown 
m Section UVDl below. 



* In Section IV A of y it is suggested that the Sultana— Dyer metric 
is equal to the McVittie metric |[T]| in which a{t) = aot^^^ and 
m(t) = mo, for some constants ao and mo (see Eq. (62) in Q). 
Let denote the latter metric by g . Indeed, since m is constant 
and in view of Proposition [T] g is conformally related to the 
Schwarzschild metric. Moreover, as one may explicitly check via 
our Eq. Illlcl l. the Einstein tensor of g has a vanishing spherical 
part. Despite sharing these two properties, g and the Sultana- 
Dyer metric are not equal. 

^ The problem with the reasoning in Section II of Q| is the following 
(numbers refer to equations in Q): It is true by construction 
that the Sultana— Dyer metric (2.1) is conformally related to the 
Schwarzschild metric, as expressed in the second line of (2.3) [the 
first line in (2.3) does not follow], but the conformal function a 
depends non-trivially on the Schwarzschild coordinates for time 
and radius (denoted by 77 and f in Q|: Cf. our discussion in 
the last paragraph of Section FlII Al l. Hence it is not possible 
to introduce a new time coordinate t that satisfies dt = adfj 
(the right hand side is not a closed 1-form), as pretended in the 
transition to (2.5). 

^ To show this, one has to show that there exists no timelike di- 
rection with respect to which the Ricci tensor (or, equivalently, 
the Einstein tensor) is spatially isotropic. This can be shown 
as follows: First note that the Einstein tensor of the Sultana- 
Dyer metric has the form Ein = fiu u + rk (g) k (see }l9|), 
where u is a normalized future-pointing spherically-symmetric 
timelike vector field and fc the in-going future-pointing light- 
like vector field orthogonal to the 50 (3)-orbits normalized such 
that g{u, k) = 1. In particular, the spherical part of the Ein- 
stein tensor vanishes: Hence, the Einstein tensor is spatially 
isotropic iff there exists a non-vanishing spacelike spherically- 
symmetric (i.e. orthogonal to the 50(3)-orbits) vector field s 
with Ein(s, s) = 0. Without loss of generality one can chose s 
to be normalized: s = sinhx"" + cosh xe, where e is the normal- 
ized vector field orthogonal to u and to the 50(3)-orbits pointing 
in positive radial direction. Hence one has k = u — e and thus: 
Ein(s, s) = sinh^(x) + T exp(2x). Clearly, the latter expression 
vanishes nowhere in the physically interesting region (cf. Eq. (26) 
in ligH ). where both fi and t are positive. 
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A. Einstein's equation for the McVittie ansatz 

In the following we will restrict to those generalizations 
of the McVittie model which keep the metric ansatz ^ 
and thus generalize only the matter model. For this 
purpose it is convenient to write down the Einstein's 
equation for an arbitrary spherically symmetric energy- 
momentum tensor T. Recall that spherical symmetry im- 
plies for the component of T with respect to the orthonor- 
mal basis (|TOl) that r(ea, e^) — and T{eA,eB) oc Sab, 
where a S {0,1} and A,Be {2,3}. Hence, the only 
independent, non- vanishing components of T are 



S 

Q 
p 
J 



T(eo,eo) 
T(ei,ei) 
T(e2,e2) 
-T(eo,ei), 



(31a) 
(31b) 
(31c) 
(31d) 



and these are functions which do not depend on the an- 
gular coordinates. Note that S is the energy density, 
Q and P the radial and spherical pressure, and J the 
energy flow — all referred to the observer field eg. The 
sing in (|31dp is chosen such that a positive J means a 
flow of energy in positive radial direction. Taking (|3ip 
into account, the Einstein equation for the McVittie 
ansatz ([Ij and an arbitrary spherically symmetric energy- 
momentum tensor T reduces to the following four equa- 
tions: 



(am)- = -47ri?2 (f)^ J 
87r5 = 3i^^ 
Svrg = -SF^ - 2^2 A 
P = Q. 



(32a) 
(32b) 
(32c) 
(32d) 



In view of ([22)l . the first equation relates the time varia- 
tion of the Weyl part of the MS energy contained in the 
sphere of radius R with the energy flow out of it. The 
last equation is nothing but spatial Ricci-isotropy. 

In the following subsections we will consider three 
models for the cosmological matter which generalize the 
original McVittie model: perfect fluid, perfect fluid plus 
heat flow, and perfect fluid plus null fluid. 



B. Perfect fluid 

Perhaps the simplest step one can take in trying to 
generalize the McVittie model is to stick to a single per- 
fect fluid for the matter, but dropping the condition Q 
of 'no-infall' by allowing for radial motions relative to 
the d/dt observer field. In this way one could hope to 
avoid a particular singular behavior in the pressure that 
may be due to the 'no-infall' condition, though it is clear 
that the persisting geometric singularity must show up 
somehow in the matter variables as already discussed in 
Section IIII Dl Unfortunately, as already shown in [3| , 
the relaxation of Q does not lead to any new solutions. 



What we want to stress here is that the reason for this, 
as shown in more detail below, lies precisely in the re- 
striction imposed by spatial Ricci-isotropy. 

We take thus the perfect-fluid energy-momentum ten- 
sor ([2]) for the matter and an arbitrary spherically sym- 
metric four-velocity u. The latter is given in terms of the 
orthonormal basis for the metric ^ by 

M = coshx eo -|- sinhx ei , (33) 

where x is the rapidity of u with respect to the observer 
held Co (a positive x corresponds here to a boost in an 
outward-pointing radial direction). The non- vanishing 
components of the matter energy-momentum tensor ^ 
with four-velocity are: 

(34a) 
(34b) 
(34c) 
(34d) 

Clearly, the case of vanishing rapidity must lead to 
the original McVittie model. In this case, in fact, the 
matter energy-momentum tensor p4p is already spatially 
isotropic so that (j32d[) is identically satisfled. Moreover, 
P2ap implies (am) =0 and hence, in view of (HH), F = 
d/a. Herewith Einstein's equation reduces to Q and 
thus one gets back the original McVittie model. 

In case of non-vanishing rapidity, spatial Ricci- 
isotropy (j32dp implies the following constraint: 



T(eo, eo) = Q+{g + p) sinh^ x 
r(eo,ei) = -(e + p)sinhxcoshx 
T(ei,ei) =p + {Q + p)fi\nh^x 
T(e2,e2) =T(e3,e3) - p. 



Q + p = Q. 



(35) 



This means that the energy momentum tensor ^ has the 
form of a cosmological constant (using ((35|) in div T = 
it implies dp = and this, in turn, using again (j35p . 
implies dg — 0) so that this case reduces to the 
Schwarzschild-de Sitter solution and hence does not pro- 
vide the physical generalization originally hoped for. 



C. Perfect fluid plus heat flow 

In a next step one may keep ([55]) and drop the con- 
dition that the fluid be perfect, in the sense of allowing 
for radial heat conduction. This is described by a spatial 
vector field q that represents the current density of heat, 
which here corresponds to the current density of energy 
in the rest frame of the fluid. Hence q is everywhere 
orthogonal to u7 The fluid's energy momentum tensor 



^ We note that the parametrization of the energy-momentum ten- 
sor given in 4] is manifestly different. Whereas we parametrized 
it in the usual fashion in terms of quantities (energy density, 
pressure, current density of heat) that refer to the fluid's rest 
system, the authors of 3 also write down I I36I I (their Eq. (79)), 
but with q orthogonal to eg (compare their Eq. (93)) rather than 
u, which affects also the definition of q. In fact, marking their 
quantities with a prime, their expression (79) is equivalent to our 
I I36II iff p = p', q = q' cosh x, and g = g' — 2q' sinh x- 
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then reads 

T = ® M + p (m ® It — g) + M (g) q + q (g) M . (36) 

Taking ([55)1 as fluid velocity and imposing the heat flow- 
vector q to be sphericaUy symmetric, we have 

q = ge g (sinhxeo + coshxei) , (37) 

where g is a function of (t,r). Note that a positive q 
corresponds to heat flowing in an outward-pointing radial 
direction. The independent non- vanishing components of 
the energy- momentum tensor are now as follows: 

r(eo, eo) = g + tanhx((£' + p) tanhx + 2(7) (38a) 

r(eo,ei) =(jf-cosh^x((£'+p)tanhx + 2g) (38b) 

T(ei,ei) = p+ isinh(2x)((^? + p)tanhx + 2q) (38c) 

T(e2,e2) =T(e3,e3) =p. (38d) 

Consider first the case of vanishing rapidity. Then the 
energy-momentum tensor is already spatially isotropic 
and Einstein's equation reduces to 

(am)' = -47ri?2g (1)2 ^gg^^ 

Sttq = 3F^ (39b) 

Sttp = -3^2 - 2F^ . (39c) 

These are three PDEs (though only time derivatives oc- 
cur) for the five functions a,m,g,p^ and q so that the 
system p9l) is clearly under-determining. However, it is 
not possible to freely specify any two of these five func- 
tions and then determine the the other three via ((39|) . For 
example, since the left-hand side of p9ap depends only 
on t, the same must hold for the r.h.s., which implies that 
q = /(i)A^(l - (TO/2r)2)2, where f{t) = -{am)- /ina^ . 
In particular, the heat flow must fall-off as l/r^. 

The easiest way to generate a solution in the case of 
zero rapidity is to specify the two functions a{t) and 
m{t), then let A,B,R,F be determined by the def- 
initions (|12I13I14[) . and finally let the Einstein equa- 
tions (I39al39bl39cp determine g, g, and p, respectively. 
Notice that if we happen to specify a and m such that 
a m is a constant, this immediately implies q = and 
F = d/a, which leads to the standard McVittie solutions. 
From (|39ap the following is evident: if g > (g < 0), that 
is for outwardly (inwardly) pointing heat flow, the Weyl 
part of the MS energy decreases (increases) , as one would 
expect. 

Now we turn to the general case with non-vanishing 
rapidity: As it was the case for the perfect fluid in the 
previous subsection, the condition (j32dp of spatial Ricci- 
isotropy implies a constraint on the matter: 

(£i + p)tanhx + 2g = 0. (40a) 

Using this, the other components of the Einstein's equa- 
tion reduces to: 

(am) ^ +4TTR^q (40b) 
ST:g = 3F^ (40c) 
&TTP = -3F2 - 2F^ . (40d) 



These are almost the same as in the case of vanishing 
rapidity (see p9|) ). except for the opposite sign on the 
right-hand side of (j40b|. This simply results from the 
fact that, according to (j38b[) . J = --T(eo,ei) = q for 
vanishing rapidity, whereas, due to the constraint (|40ap . 
J — — r(eo,ei) = — g for non- vanishing rapidity. This 
will be further interpreted below. Notice that for the 
equation of state g+p = (cosmological term) (|40ap im- 
plies g = 0, thus leading once more to the Schwarzschild- 
de Sitter solution (see comment below Eq. ([35]) ). Hence- 
forth we assume g+p ^ 0, which implies that one can 
solve the constraint (|40ap for the rapidity: 

tanhx = ^ (41) 

g + p 

provided that \2q/{g + p)\ < 1. 

The Einstein equation gives now four equations for the 
six functions a,m, g,p,q, and x- As in the case of van- 
ishing rapidity, this system is under-determining and it 
is not possible to freely specify any two of these six func- 
tions and then determine the the other four. In a similar 
fashion as before, the easiest way to generate a solution 
is to specify the two functions a{t) and m{t), to let then 
the definitions (|12ll3ll4p determ ine A,B,R,F , and fi- 
nally use the Einstein equations (|40bl40cl40dp and (|4T|) 
to determine q,g,p^ and Xi respectively. Again, choos- 
ing a and m such that their product is constant implies 
g = and F ~ a /a, which leads to the standard McVittie 
solutions. 

In passing we remark that the condition g + p > can 
be expressed geometrically in terms of the second time- 
derivative of the areal radius. Indeed, adding either (|39bp 
to (I39cl) or gOcl) to jiOd)) we obtain, taking into account 
eo = (A/B)d/dt and US]): 

MP+p) = -e„(^), (42) 

which is positive iff the rate of change Cq [R) /R is a de- 
creasing function along the integral lines of the observer 
cq. In other words, g + p is positive iff ln(i?) is a con- 
cave function on the worldline of the observer bq, which 
is implied by, but not equivalent to, the function R being 
concave. 

From (|4ip and (j40bp . and assuming g + p > 0, one sees 
the following: If x > (x < 0), that is for an outwardly 
(inwardly) moving fluid with respect to eg, we have g < 
(g > 0), that is an inwardly (outwardly) pointing heat 
flow, and the Weyl part of the MS energy decreases (in- 
creases). This means that the heat flow's contribution 
to the change of E„ never compensates that of the fluid 
motion, quite in accord with naive expectation. Below 
we show that for small rapidities the contribution due to 
the heat flow is minus one-half that of the cosmological 
matter. 

Let us now return to the sign-difference of the right- 
hand sides of (|39ap and (j40bp . From ([55]) one infers that 
J is the sum of the two contributions coming from the 
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heat flow 

(7(1 + 2 sinh^x) (43) 

and from cosmological matter 

Jm {g + p)smhxcoshx, (44) 

respectively. The constraint (|40ap can be written in the 
form 

2cosh2x-^h + (l + 2sinh2x)-^m = 0, (45) 

which, for smaU rapidities x (that is neglecting quadratic 
terms in x), impUes 2Jh + Jm ~ 0. In this approximation 
the spatial energy-momentum flow due to heat is minus 
one-half that due to the cosmological matter. For the 
total flow this implies J = Jm + Jh ~ Jm/'i ~ —Jh- Now 
the sign difference between p9ap and (|40b[) is understood 
as follows: In case of vanishing rapidity one has Jm — 0, 
Jfi = q and hence J — q (leading to (|39ap ). whereas 
a short calculation reveals that in case of non-vanishing 
rapidity the constraint (|45|) implies J — Jm + Jh = —q, 
leading thus to (|40bp . 



D. Perfect fluid plus null fluid 

The last tentative generalization we consider is tak- 
ing for matter the incoherent sum (meaning that the 
respective energy-momentum tensors adds) of a perfect 
fluid (possibly with non-vanishing pressure) and a null 
fluid (eventually representing electromagnetic radiation) . 
This clearly contains as special case the matter model 
considered by Sultana and Dyer [l9l | in which the pres- 
sure vanishes. We already stressed in Section fill Al that 
the metric ansatz of fl^ is different from ([1]). Here we 
show that the matter model of 19] is essentially incom- 
patible with ^ except for trivial or exotic cases. 

The matter model consists of an ordinary perfect fluid 
and a null fluid (e.g. electromagnetic radiation) with- 
out mutual interaction. Hence the matter's energy- 
momentum tensor is just the sum of ([2]) and 

r± = A2i±®i±, (46) 

where A is some non-negative function of t and r and 
Z+ and l~ are, respectively, the outgoing and ingo- 
ing future-pointing null vector fields orthogonal to the 
spheres of constant radius r partially normalized such 
that g{l~^,l~) = 1. (It remains a freedom i-^ a^^l^, 
where a is a positive function) . Without loss of generality 
we make use of this freedom and choose: 

Z± = (eo±ei)/V2, (47) 

where Bq and ei are the vectors of the orthonormal 
frame (|10|) . The components of the whole energy- 
momentum tensor with respect to this frame are then: 



T(eo, eo) = g+{g + p) sinh^ X + 5 (48a) 

T(eo , ei ) ^ -{g + p) sinh x cosh X T ^ (48b) 

T(ei,ei) =p+ (i? + p)sinh2x+ (48c) 

r(e2,e2) =T(e3,e3) (48d) 

Here and below the upper (lower) sign corresponds to the 
outgoing (ingoing) null field. 

In the present case, the condition (j32dp of spatial Ricci- 
isotropy is equivalent to the constraint: 

{g + p)smh\ + ^X^ ^Q. (49) 

In the physically relevant case in which g + p > 
this equation has only the trivial solution x = and 
A = 0, which leads to the original McVittic model. In 
the case g + p — (l49|l implies A = 0, leading thus 
to the Schwarzschild-de Sitter spacetime (see comment 
below (|35p). Hence, a new solution is only possible if 
the matter is of an exotic type that satisfies g + p < 0, 
which either violates the weak energy-condition (g > 0), 
or, less catastrophically, the dominant-energy condition 
(g > \p\). In particular, for the matter model considered 
by Sultana and Dyer, one would need to violate the weak 
energy-condition. 

V. CONCLUSION 

We conclude by commenting on the the main dif- 
ferences between these generalizations and the original 
McVittie model. First we stress once more that neither 
allowing for a nonzero rapidity nor a nonzero heat flow 
can eliminate the sing:ularity at r — ni/ 2 {R — 2am) (as 
erroneously stated in p!|). The only substantial new fea- 
ture of these generalizations is that the Weyl part of the 
MS energy E„ — am is not constant anymore. In view 
of the fact that the combination 

m/r = A^E^/R « E„/R (50) 

contained in the McVittie ansatz gives the 'Newtonian' 
part of the potential in the slow-motion and weak-field 
approximation (see [H), we deduce that in order to get 
the geodesic equation for the generalized McVittie model, 
it suffices to substitute mo with in the equation of 
motion derived in |l| . This means that the strength of the 
central attraction varies in time according to (|32ap . lead- 
ing to an in- or out-spiraling of the orbits if di?w(eo) > 
or dE„{eo) < 0, respectively. 

We identified the origin of why we could not vary the 
rapidity and the heat fiow independently in the condi- 
tion (|40ap of spatial Ricci-isotropy, which is built into 
the ansatz ([1]). We saw that this geometric feature ren- 
ders this ansatz special, so that it would be improper to 
call it a general ansatz for spherical inhomogeneities in a 
fiat FLRW universe. It remains to be seen whether use- 
ful generalizations exist which are captured by equally 
simple ansatze. 
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APPENDIX A: PROOF OF PROPOSITION [1] 

In this appendix we compute the intersection of the 
set 5mcv of metrics of type ([1]), which we denote in the 
following by g"'^, with the set iS^s of metrics conformally 
related to an exterior Schwarzschild metric. Explicitly, 
the latter are of the form 

(Ala) 

where 

denotes the Schwarzschild metric with mass Mq in 'stan- 
dard' coordinates. The question is: for which functions a 
and m and, respectively, for which function and pa- 
rameter Mo does the equation g^"^ = ^"^dtt™ hold? 
Such an equality can be eventually established by finding 
a coordinate transformation, cj) say, between the coordi- 
nates^ [t, r) in H]) and (T, R) in (|XT|) which brings ^ 
in form (jAip . This involves solving coupled, non-linear 
partial differential equations for 0, which depend on the 
four unknown parameter a,TO, il, and Mq. Needless to 
say that this is not really a thankful task. Alternatively, 
a better approach would be to compare all the indepen- 
dent, algebraic curvature-invariants of the two metrics: 
This would lead to a system of equations between scalars 
which involves the coordinate transformation in an al- 
gebraic way (i.e. non differentiated). 

We adopt here an approach which is somewhere in the 
middle: First, we use just three invariants (the areal ra- 
dius and the Ricci and the Weyl part of the MS energy) 
to drastically restrict the form of the coordinate transfor- 
mation (see (jA14l) ) and derive thereby constraints on the 
free parameters a, m, 51, and Mq (see (jA10[) and (jA13l) ). 
Second, we perform this restricted coordinate transfor- 
mation and determine it completely. To simplify the 
calculation, instead of g"!"^ — H.'^g^^^, we consider the 
equivalent equation Q^^g^^^^ = dlj™- In fact, for the 
Schwarzschild metric (jAlb[) it is immediate that the 
above mentioned quantities are, respectively, given by: 

RiolT) (A2) 
E^{gl:^-)^Mo, (A3) 
i^a«")=0. (A4) 



The transformation between the angular variables is just the 
identity. 



In order to compute the respective quantities for the met- 
ric ri~^g"^ we first give their scaling behavior under 
conformal transformations. 

Clearly, because of their very definitions, for the areal 
radius and the Weyl part of the MS energy it holds: 

Ri^^g) = *i?(sf) (A5) 

and 

£;w(*'g) = *£^w(g), (A6) 

respectively. For the whole MS energy it easily follows 
from (HOI) and ([Mjl : 

E{^^g) = * (E(g) + R^g{VR, Vln *) 
-Fii?3£/(Vln«', Vln*)j , 

where all the quantities on the r.h.s. are referred to the 
metric g. Hence, taking the difference between (IA7p 
and (jA6|) one gets that the Ricci part of the MS energy 
scales exactly like the whole MS energy, that is according 
to (lAT)) . 

Using these scaling properties together with (fT3|) 
and (|22p we get immediately: 

m^'OaZ) =n-\l + m/2rfar (A8) 
E^in~'g:%) = n-'am. (A9) 

The equality between the Weyl part of the MS en- 
ergy (K9\i and (K3\i implies 

fl{t, r) Mo = a{t) m{t) , (AlO) 

which gives a condition between the parameter a, to, f2, 
and Mq. Since we assumed that Mq is positive, ()A10|) 
can be read as the expression for the conformal factor in 
the {t, r) coordinates. This, together with the equality 
between the areal radius (jASp and (|A2p . implies in turn 

R{t,r) = ^{\ + mit)/2rfr, (All) 
7n(t) 

which gives the first component of the coordinate trans- 
formation (j). Now, using the scaling property (|A7p for the 
Ricci part of the MS energy, the expressions (|2ip and (fT5)) 
for the Ricci part of the MS energy and, respectively, the 
areal radius of the McVittie metric ansatz, and (|A10p for 
the conformal factor, one gets, after some computations, 

E^in-^g^Z) ^ ^{A^arf . (A12) 

2am \ / 

The equality between (jA12p and (|A4p then implies 

TO = , (A13) 

that is TO = Too for some positive constant toq. This, in 
turns, implies that the transformation (jAlip for R de- 
pends only on r and not on t. Since the metrics are both 
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in diagonal form, this implies that the transformation for 
T must depend on t only. 

Summarizing, so far we have seen that a set of neces- 
sary conditions for the equality of the two metrics implies 
the constraints ()A10P and ()A13|) and that the coordinate 
transformation between (t, r) and (T, R) is of the form 



m = fit) 



Rir) 



Mo 



{l + mo/2rfr, 



(Al4a) 
(A14b) 



for some differentiable function / of t. Now, explicitly 
expressing the metric fi^gr^j^" in the {t, r) coordinates 
according to the coordinate transformation (|A14p and 
the constraints (jAlOp and (|A13p . and putting the result 
equal to g"^, the only new condition that one gets is 



mo a 



(A15) 



Here, the plus can be chosen in order to exclude a time 
inversion. It is important to note that (jA15|) (together 
with an initial value) determines / uniquely and do not 
give any constraint on the parameters a,m, f2, and Mq: 
The only constraints remain thus (jAlOp and (jA13p . 

The proof is concluded noticing that (|A10[) means that 
the only constraint on is that, expressed in the {t,r) 
coordinates, it depends on t only and hence, in view 
of (|A14a|) and expressed in the (T, R) coordinates, that 
it depends on T only. More geometrically, this can be 
restated saying that the gradient of Q must be propor- 
tional to d/dT, the Killing field of the Schwarzschild 
metric (see (jAlbp ). 



APPENDIX B: PROOF OF PROPOSITION [2] 



Inserting the definition ([14]) of F in the expression ([231) 
for the Ricci scalar and organizing the result in powers 
of rB = r — ni/2 we get: 



Seal = - 12 I - 
a 



6 / , a(a m) ,1a 
rB \ 



(rBy 



[a m) 



(am) a{am) 



{rB) 



3 ,m(am) 
rA — - 



(Bl) 



Hence, the Ricci scalar remains finite in the limit r 
m/2 iff all the three coefficient of {rB)~'', for k e 



{1, 2, 3}, vanish in this limit, that is iff it holds: 

d(a to) 



(a to) 
a 

m(a to) — . 



r 


)- 


TO — 






2 


' (a to) 


+ TO 




a 



(B2a) 



= (B2b) 
(B2c) 



These conditions are clearly understood to hold for all 
times t in which the functions a and to and their deriva- 
tive exist. In view of ([B2cP we have to distinguish be- 
tween two cases: to = and (a to)' = 0, respectively. In 
the first case the system ([B2p reduces to the set of con- 
ditions m{a/a + 3(d/a)^) — and m{d/a+ (d/a)^) = 0, 
which, in turn, reduces to to = (and a arbitrary), corre- 
sponding to the FLRW metric, or to d = (and to = 0), 
corresponding to the Schwarzschild metric. In the second 
case, in which (a to)' = 0, ([B2p reduces to m{a/a)' — 0, 
which implies either to = (and a arbitrary) , correspond- 
ing again to the FLRW metric, or {a /a)' = 0. Together 
with (a to)' = 0, the latter corresponds to a McVittie 
metric with exponentially-growing (or -falling) scale fac- 
tor a{t), that is to a Schwarzschild-de Sitter metric. 



APPENDIX C: SHEAR-FREE OBSERVER 
FIELDS IN SPHERICALLY SYMMETRIC 
SPACETIMES 

Towards the end of Section[ni we made use of the fol- 
lowing result: A spherically symmetric normalized time- 
like vector field m in a spherically symmetric spacetime 
(Al , g) is shear free iff the metric hu that g induces on 
the subbundle u-^ := {v e TM \ g{v, u) — 0} by restric- 
tion is conformally flat. 

To prove this, we first note that the subbundle is 
integrable, in other words, u is hypersurface orthogonal. 
This follows from the spherical symmetry of u, which 
implies that contains the vectors tangent to the 2- 
dimensional 50(3) orbits. Hence u essentially lives in 
the 2-dimensional orbit space^, where it is trivially hy- 
persurface orthogonal. The hypersurfaces orthogonal to 
u in 4-dimensional spacetime are then the preimages un- 
der the natural projection of the hypersurfaces (curves) 
in the 2-dimensional orbit space. 

As a result, we may now locally introduce so-called 
isochronous comoving coordinates, with respect to which 
u = A{t,r)-^d/dt and 



g^A^{t,r) dt" -B\t, r) dr^ - R^t, r) g^ 



(CI) 



^ The orbit space is the quotient M/ ~, where ~ is the equivalence 
relation whose equivalence classes are the orbits. It is a manifold 
on the subset corresponding to 2-sphere orbits, to which we re- 
strict attention here. To say that "it essentially lives in the orbit 
space" means that it is the pull-back of a 1-form on the quotient 
via the natural projection. 
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(Note the different meanings of the functions A and B as 
compared to (HH)). We now consider the tangent-space 
endomorphisms Vm : X Vx^ and their projection 
into the orthogonal complement of u, i.e., 

V^u:={P^oVuoP^)\^^, (C2) 

where id — u (E) u is the projection orthogonal 

to u (id is the identity endomorphism in the tangent 
spaces of A^). Note that V^tt is symmetric due to the 
hypersurface orthogonality of m. A direct computation 
using (jCip yields 

V^u = u{ln{B)) Pr + w(ln(i?)) P52 , (C3) 

where P,. and P52 are the projections parallel to d/dr 
and parallel to the tangent 2-planes to the S'^-orbits, re- 
spectively. The trace 9 of V^m, which gives the expan- 
sion of M, is 6* = u{\n{B)) +2M(ln(i?)), so that the trace- 
free part of V^u, known as the shear endomorphism cr, 
is given by: 

(T := V^M - ie* id-L = a {Ps2 - 2Pr) , (C4) 

where a := ^u(ln.{R/B)^ denotes the shear scalar (only 
defined in a spherically-symmetric setting) and id^ = 
Pr + Ps^ the identity endomorphism in u^. In passing, 
we note that the defining equations for 9 and a just given 
immediately lead to the following simple relation between 
the shear scalar, expansion, and the variation of the areal 



radius along u, that we made use of in Section lTlI Bl 

cr + 0/'3 = u(ln(i?)) . (C5) 

Now, according to (|C4[) . the shear of u vanishes iff the 
shear scalar a does, that is, iff u(ln{R/B)) vanishes. 
This is equivalent to R/B being independent of t or to 
R{t,r) = iJ.{r)B{t,r) for some function /i, so that the 
line element (jCl|) can be rewritten in the spatially con- 
formally flat form 

g^A^it,p) dt" -C^{t, p) {dp" + p2 g^, ) , (C6a) 
where A{t,p) := A{t,r{p)), C{t,p) := C(t, r(p)), and 

Cit,r) = ^Ml^ (C6b) 

with 

pir)=poe.^[l^^]. (C6c) 

Hence we see that vanishing shear of u implies confor- 
mal flatness of the corresponding spatial metric. For the 
converse we first note that, since u and g are spherically 
symmetric, the spatial metric is itself spherically sym- 
metric, so that g can be written in the form (jC6ap . This 
implies that the corresponding R/B depends only on the 
radial coordinate and hence that the shear of u vanishes. 
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